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Abstract 

We study a multispecies one-dimensional Calogero model with two- and three-body 
interactions. Using an algebraic approach (Fock space analysis), we construct lad- 
der operators and find infinitely many, but not all, exact eigenstates of the model 
Hamiltonian. Besides the ground state energy, we deduce energies of the excited 
states. It turns out that the spectrum is linear in quantum numbers and that the 
higher-energy levels are degenerate. The dynamical symmetry responsible for degen- 



eracy is SU(2). We also find the universal critical point at which the model exhibits 
singular behaviour. Finally, we make contact with some special cases mentioned in 
the literature. 
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1 Introduction 



The ordinary Calogero [1] model describes N indistinguishable particles on the line 
which interact through an inverse-square two-body interaction and are subjected to 
a common confining harmonic force. The model is completely integrable in both the 
classical and quantum case [2], the spectrum is known and the wave functions are 
given implicitly. The model and its various descendants (also known as Calogero- 
Sutherland-Moser systems [3]) are connected with a number of physical problems, 
ranging from condensed matter physics [4] to gravity and black hole physics [5] . The 
algebraic structure of the Calogero model, studied earlier using group theoretical 
methods [2,6], has recently been reconsidered by a number of authors in the frame- 
work of the Sn (permutational) algebra [7]. This operator approach is considerably 
simpler than the original one, yields an explicit expression for the wavefunctions 
and emphasizes the interpretation in terms of generalized statistics [8], especially 
Haldane's exclusion statistics [9]. In Haldane's formulation there is the possibility 
of having particles of different species with a mutual statistical coupling parameter 
depending on the i th and j th species coupled. On the level of the Calogero model, 
this corresponds to the generalization of the ordinary ID Calogero model with iden- 
tical particles to the ID Calogero model with non-identical particles. This can be 
done by allowing particles to have different masses and different couplings to each 
other. In this way we obtain a ID multispecies Calogero model. Very little is known 
about its spectra and wavefunctions [10-13]. 

In the present Letter, which is in a sense a continuation of our investigation of 
the ordinary Calogero model [14], we use an algebraic (operator) method to find 
some of the salient features of the multispecies Calogero model on the line with two- 
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and three-body interactions. In Section 2 we present the Sn invariant Hamiltonian 
H of the model, together with its ground state wavefunction and ground state 
energy. After performing a certain similarity transformation, we get a much simpler 
Hamiltonian H, which we separate into parts describing the center-of-mass motion 
and the relative motion of particles. We express H in terms of generators of the 
SU(1, 1) algebra. All analysis is made in Hilbert space. Section 3 contains our most 
important results. By applying Fock space analysis, we find some of the excited 
states of H, their energies and degeneracies. Closer inspection of the Fock space 
that corresponds to the relative motion of particles reveals the existence of the 
universal critical point at which the system exhibits singular behaviour. This result 
generalizes that mentioned in [14] . We also establish the necessary conditions for the 
equivalence of the two multispecies Calogero models. In Section 4 we briefly repeat 
the main points of the paper and make contact with the models studied in [10-13]. 
We particularly discuss the necessary conditions for vanishing of the three-body 
interaction in the starting Hamiltonian H. 

2 A multispecies Calogero model with a three- 
body interaction 

Let us consider the most general Calogero type ground state for the iV-body quantum 
mechanical problem on the line (h — 1): 

^o(x 1 ,x 2 , ■■■x N )= Ae-?^i^, (1) 

where the prefactor A is given by 

A = - x 3 ) Vii , Vij = vji, i,j — 1,2 - ■• N 

i<j 



3 



and Uij are symmetric statistical parameters between particles The harmonic 

frequency uj is taken to be the same for all particles. Masses of the particles ( rrii ) 
are, in general, not equal. 



The Hamiltonian which possesses the above state (1) as the ground state is 
a 

dxi 



^ i,j,k£ m j( x j ~~ x i)( x j ~ x k) 



HV (x 1 ,x 2 , ■ ••xn) = E ^o(x 1 ,x 2 , ■ ■•xn), 

E = uj(— + -2^v ij ). 

The symbol (i,j,k ^) in the last term denotes the summation over all triples of 
mutually distinct indices. 

The Hamiltonian in Eq.(2) is invariant under the group of permutation of N ele- 
ments, Sn- The elementary generators Kij of the symmetry group Sn exchange 
labels i and j in all quantities, according to the rules: 

KijXj KijTflj TTij-Kjj, KijVji Vi\Kij, 

Kij Kji, [Kij} 1, 
KijKji = KjiK a = KuKij, for i^j,i^l,j^ I. 



It follows that KijHKij = H, i.e. [H, K i3 \ = 0, Vi, j. 

In a sense, one may think of (2) as a smooth, sufficiently small deformation of the 
ordinary Calogero model. 

A few additional remarks concerning the Hamiltonian (2) are in order. 

(1) It describes distinguishable particles on the line, interacting with harmonic, two- 
body and three-body potentials. As far as we know, it is the first time that such 
Hamiltonian is considered in full generality. An earlier attempt to solve the similar, 
but less general Hamiltonian ( with all masses rrii equal) can be find in [11] (see also 
Ref.[13]). 

(2) The asymptotic behaviour of its eigenstates should be ^ oc (x { — Xj) Uij as 
(xi - Xj) ->• 0. 

(3) A well-known stability condition demands that the two-body couplings 
Vijiyij — 1) should be greater than — |, Vi,j. 

(4) Setting i/^ = u, and rrii = m, Vi, we recover the ordinary N-body Calogero 
model [1]. In that case, owing to the identity J2cyci. ( x --x-)(x--x k ) = ® which holds in 
ID, the three-body term in Eq.(2) trivially vanishes. When are mutually different 
but rrii = tti, Vi, we recover the model treated in [11,13]. Finally, when i/^ = arriimj, 
a being constant, we obtain the model mentioned in [10,12]. We comment on these 
cases in Section 4. 

In the following we analyse the most general case, namely the Hamiltonian of Eq.(2) 
without restrictions (4), given above. 

Let us perform the non- unitary transformation on ^o(xi, X2, ■ ■ • Xn)'- 

V (x 1 ,x 2 ,---x n ) = A- 1 V (x 1 ,x 2r --x N ) = e-i^i^. (3) 

It generates a similarity transformation which leads to another Sn invariant (but 
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non-Hermitean) Hamiltonian H: 

H = A~ l HA, 

H^ (x 1 ,x 2 , ■■■x N ) = E ^ (x 1 ,x 2 , ■ --xn), 

with Eq given in Eq.(2). 
We find H as 

2 rrn dx\ + 2 ^ miXt 2 ^ (a* - xj) ( m, ftr, m, dx 3 h U 

Notice that in Eq.(4) two-body and three-body interactions apparently disappeared 
but they are hidden in the last term of Eq. (4). 
It is convenient to introduce the variables (X, £j) 

^= y^ ttE^ > 6 = ^-X, < = 1,2,...JV, (5) 

2^i=l m i M i=l 

and the linear combinations of derivatives (^-, ^) 

A JL-H!iJL l = f 1 - = 19 /V 

ae, ar, Max ' ax ^a**' 1 i ' / '"' iV - 

Note that the variables & ( as well as ) are not linearly independent, i.e. 

In terms of the variables just introduced, the Hamiltonian H, Eq.(4), separates into 
parts which describe its center-of-mass motion (CM) and its relative motion (R), 
namely H = H C m + H R , with 

Id 2 1 9 P 

1^1 (9 2 1 ^ 1 ZA 1 (9 1(9 
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The Hamiltonian H, Eq.(2), can also be decomposed into Hqm and Hr. 
The wave function (3) separates as 

*o(*i, x 2 , ■ ■ ■ x N ) = (X)M^, 6 • • • Zn) = e~^ x2 e-% ^ "< 

The ground state energy E splits into the energy of CM ( E Q cm = \w ) and the 
energy of relative motion (E 0R = ^f^-u> + \ujY,i^j v ij)- 
We define the set of operators {T + ,T_,T } as 

1 N 

T + =oJ2 m i X l 
Z i=l 

T _ 1 y 1 92 | 1 y ( 1 g 1 g ) 

2^ 1 m i dx 2 2f£ j (x i - Xj) rriidxi rrij dxj 

^(g^ + f)- (8) 
Using Eqs.(5) and (6) one can easily show that these operators also split as 
T±,o = T±,o (CM) + T±,o (R) ■ 

Operators {T + ,T„,T } satisfy the SU (1,1) algebra: 

[T_,T + ]=2T , [T ,T ± }=±T±. 
In terms of these operators, the Hamiltonian (4) reads H = u 2 T + — T_. 

3 Ladder operators and Fock space representa- 
tion 

Now we introduce pairs of creation and annihilation operators {Af, A±} and {A 2 , A 2 }: 



Af = l(uT + + -T^TT , (9) 
which satisfy the following commutation relations: 

[AZ,A+] = 1, [A^,At] = -H, 

CO 

[■A-i , A 2 ] = [A*, A%\ = 0, [A-y , A£] = Af, [Af, A 2 ] = —A x , 

[H,Af]=±uAf, (10) 
[H,A£]=±2uAf. 
They act on the Fock vacuum |0) oc ^> (xi,x 2 , • ' ' x n) as 

a 1 -\o) = a 2 \o)=o, (o|o) = i. 

The excited states are built as 

\n u n 2 ) oc A+"M+" 2 |0), Vm, n 2 = 0, 1, • • • (11) 

The repeated action of the operators Af ( A 2 ) on the vacuum |0) reproduces 
in the coordinate representation Hermite polynomials ( hyp ergeo metric function ), 
respectively. 

The states \rii,n 2 ) are eigenstates of H, Eq.(4), with the eigenvalues 
( cf. Eq.(10) ) 

E nuna = u(m + 2n 2 ) + E . (12) 

The energy spectrum is linear in quantum numbers ni, n 2 . This result is universal, 
i.e. it holds for all parameters and Vij in the Hamiltonian H (or H ). Notice 
that the energy of the ground state and the energy of the first excited state are 
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non-degenerate whereas the higher energy levels are degenerate. The structure of 
degeneracy is as follows: 



n-y 
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n = n\ + 2ri2 


Degenerate states 











|0> 


1 





1 




2 





1 


2 
2 


Af|0) 


1 

3 


1 




3 
3 


a+a+|o> 

A+ 3 |0) 



2 
4 


2 
1 



4 
4 
4 


A+ 2 A+|6> 
<|0) 


5 
3 
1 



1 
2 


5 
5 
5 


Ario> 

At 3 At\0) 
AfAt 2 \0) 











It is evident that for n = even, the degeneracy is (| + 1) and for n = odd, the 
degeneracy is ( z? y i ). 

In order that the two models described by Hamiltonian (2), with statistical param- 
eters vij and v^, have the same tower of states (11) and the same spectrum (12), 
the necessary conditions are J2i<j v %j — Si<j v 'iji u = oj' and the number of particles 
should be the same, i.e. iV = N'. 

The dynamical symmetry algebra of the model is defined as maximal algebra com- 
muting with the Hamiltonian H. The dynamical symmetry of the ordinary Calogero 
model is complicated polynomial algebra denoted by CV(^) in [15]. In our case, ow- 
ing to the fact that H (10) can be rewritten in terms of two independent, uncoupled 
oscillators (see bellow Eqs.(13) and (14)), this polynomial algebra can be linearized 
to the ordinary SU (2) algebra. This is the minimal symmetry that remains in the 
generic case, i.e. for general Uij and rrii. In fact, this is the same dynamical symmetry 
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underlying the two-body Calogero model [15,16]. 

We point out that one can construct the creation and annihilation operators 
W,5 2 -}: 

Bf = A 2 ± - \A?. (13) 

In terms of Af and Bf, the above-mentioned SU{2) algebra ([J + , J_] = 2J , 
[J , J ± ] = ± J ± ) is generated by 

1 



J+ = A+ 2 5 2 



^2(^-1 + ^)^ + 1) 
1 



J = ^(N 1 -2N 2 ). 

Here, iVi and N 2 are number operators counting A\- and £> 2 - modes, respectively. 

The benefit of the construction (13) is that the operators Af, corresponding to 
the CM motion, decouple completely (cf. Eq.(lO)), i.e. 

[Af,B?\=0. 



Hence, we get 



H C m = tM A i~, A t}+ = ^( N i + 



2 v x ^ 

# i? = a;[ J B 2 -, J B 2 + ] =cj(2iV 2 + 



a; 

?±1 LO, , D± 



[^,^ = ±2^. (14) 
The Fock space now splits into the CM-Fock space, spanned by Af ni \0) C M, and the 



R-Fock space, spanned by S 2 \0)r, where \0)cm oc ^o(X) and |0)/j oc ^0(6)^2 • • •£ 
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At this point it is useful to make a contact with the Fock space of the ordinary 
Calogero model. In the ordinary Calogero model the SAT-symmetric subspace of 
the whole Fock space is spanned by the states ,4+ ni |0)cM and (B 2 n2 ■ ••Bp lN )\0) R 
[7,16], where 



N N i 

A = E Oi, &k = - ^?A) fc , 2 < k < N. 
i=i i=i JV 

The so-called Dunkl-Polychronakos operators a« satisfy the algebra [7] 

N 

[cn, a]} = 8ij(\ + v E - f-Kij, (15) 
fc=i 

[a},ot] = [a^dj] = 0. 

Our operators A 1 and i? 2 correspond exactly to the operators Ai and B 2 in the or- 
dinary Calogero model. Within our algebraic treatment, we are unable to construct 
the eigenstates of (2) which correspond to the Calogero-states (B^ U3 ■ ■ ■ B^f nN )\0) R . 

Reducing the problem (2) to the (4) and (14), i.e. H — > H — > H R , has an 
interesting consequence, namely the existence of the universal critical point, defined 
by the null-vector 

umB}\0U = N-l 1 

r(0\0)r 2 2^" y ' 

The above relation (16) follows directly from Eq.(14) by demanding 
r(0\Hr\0)r = E 0R = 0. More generally, from Eq.(14) (see also Ref.[17] ) immedi- 
ately follows that 

B+B 2 = N 2 (N 2 - 1 + ^) = 0(A> 2 ), 



UJ 



B 2 B+ = (N 2 + 1)(A> 2 + — ) ee 0(iV 2 + 1), 

UJ 



*{o\BrB$-\o) R _fi m (17) 



r(0\0)r 

11 



Since Eq.(16) implies that 0(1) = ^ a = 0, it is obvious that the critical point 
(Eq.(16)) is unique, i.e. there are no similar critical points when norms of states 
involving higher powers of the operators B 2 are involved. 

At the critical point the system described by Hr collapses completely. This 
means that the relative coordinates, the relative momenta and the relative energy 
are all zero at this critical point. There survives only one oscillator, describing the 
motion of the centre-of-mass. This singular behaviour was first noticed in [14] for 
the case i/^ = v and mi = m. Of course, for the initial Hamiltonian H, which is 
not unitary ( i.e. physically ) equivalent to H, this corresponds to some < and 
the norm of the wave function (1) blows up at the critical point. For negative 
but greater than the critical values (16), the wave function is singular at coincidence 
points but still quadratically integrable. 

4 Discussion and outlook 

In this Letter we have studied the most general multispecies Calogero model on 
the line, Eq.(2) , with a three-body interaction and an extended SV invariance. 
By applying the similarity transformation A, we have obtained the Hamiltonian 
(4), on which we have performed the Fock space analysis (9,10) and found some 
of its (but not all) excited states, Eq.(ll), and their energies E nim , (12). It turns 
out that the energy (12) is linear in quantum numbers n\ and n 2 and there is a 
dynamical SU(2) symmetry responsible for the degeneracy of higher-energy levels 
with n = rii + 2n 2 > 2. By splitting the Fock space into the CM-Fock space and 
the R-Fock space (14), we have detected the universal critical point (16) at which 
the system exhibits singular behaviour. 
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To conclude this analysis, let us consider the last term in (2) more closely, namely 



^ i,j,k^ m '3 ( X 3 ~ X i){ X j ~ X k) 

If we put nij = m = const, in (18), Vj, symmetrize under the cyclic exchange of 
the indices (i — > j — > k — > i) and reduce the sum to a common denominator using 
the identity 

Et— T=0. 



cycl. 



X 



x j){, x i x k) 



we obtain that the necessary condition for vanishing of the three-body interaction 
is Vij = v = const., In this way, the problem (2) is reduced to the ordinary 

N-body Calogero model with two-body interactions only [1]. 

For the general Vij and mj, the above procedure yields the following necessary con- 
ditions for the absence of the three-body interaction (18): 



L'ijfjk ^jk^ki ^ki^ij 



V(i,j,k) (19) 



nij nik nii 

The unique solution of these conditions is Vij = a mi mj, a being some universal 
constant. This particular connection between masses and interaction parameters 
was also displayed in [10,12]. In [12], the condition (19) arose from the demand 
that the asymptotic Bethe ansatz should be applicable to the ground state of a 
multispecies many-body quantum system obeying mutual statistics, while in [10] its 
origin was not obvious. In our approach, it has the simplest possible interpretation. 
The results presented here are easily applied to the model with F distinct families 
of particles. For example, for a — 1, the ground state energy becomes 

N * 2 N a (N a -l) 1^ 
= ^(t- + 2^ m a 5 + « 1^ rn a m b N a N b ), 

A a=l A A a^b 
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F 

N =J2 N a- (20) 

a=l 

Two systems characterized by {u, m a , N a } and {uj' , m' a , N' a } are identical if X) N a = 
EK,u = uj' and = K- 

The open problem that still remains is the construction of (generalized) Dunkl- 
Polychronakos operators and a\ (15), which may help in finding the complete set 
of eigenstates of the Hamiltonian (2) (or (4)), in one-to-one correspondence with the 
ordinary Calogero model. 
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